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We study the zero-temperature transport properties of one-dimensional normal metal- 
superconductor (NS) junctions with topological superconductors across their topological transi¬ 
tions. Working within the Blonder-Tinkham-Klapwijk (BTK) formalism generalized for topological 
NS junctions, we analytically calculate the differential conductance for tunneling into two models of 
a topological superconductor: a spinless intrinsic p-wave superconductor and a spin-orbit-coupled 
s-wave superconductor in a Zeeman field. In both cases we verify that the zero-bias conductance 
is robustly quantized at 2 e 2 /h in the topological regime, while it takes nonuniversal values in the 
nontopological phase. The conductance spectra in the topological state develops a peak at zero bias 
for certain parameter regimes, with the peak width controlled by the strength of spin-orbit coupling 
and barrier transparency. 
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I. INTRODUCTION 

The search for topological superconductors in solid- 
state systems is motivated by the possibility of realizing 
Majorana zero-energy modes at their surfaces, which are 
of both fundamental and technological importance 1 ®. 
In the absence of intrinsic topological superconductors, 
much effort has been put into engineering such systems 
from conventional components 4 ®. In addition to the 
proposals involving semiconductor-superconductor hy¬ 
brid structures'® as hosts for Majorana modes, which 
have attracted considerable experimental attentiorP®®, 
there have been many recent theoretical proposals for ar¬ 
tificially engineering effectively spinless low-dimensional 
p-wave topological superconductorJ^Ml, which could lo¬ 
calize zero-energy Majorana modes at suitable defects 
such as vortex cores or system boundaries. The sub¬ 
ject has been extensively reviewed in the recent litera- 

turepnHS) 

It is particularly desirable to realize spiuless p-wave su¬ 
perconductors, as they support a single Majorana mode 
at their boundaries 1 2 . The most promising of these pro¬ 
posals involves proximity-inducing superconductivity in a 
spin-orbit-coupled semiconducting nanowire in the pres¬ 
ence of a magnetic field®®, which has subsequently been 
the subject of a number of experiments 9 " 15 . By vary¬ 
ing the magnetic field, the system is predicted to un¬ 
dergo a transition from a nontopological to a topological 
phase. Such an external magnetic-field-induced topologi¬ 
cal quantum phase transition has the considerable advan¬ 
tage of tuning the existence (or absence) of the Majorana 
mode in the experimental system simply by changing the 
Zeeman field. A key signature of the topologically non¬ 
trivial state is the quantized value 2e 2 //i of the differen¬ 
tial conductance for tunneling into the wire at zero-bias 
voltage. This quantized conductance, associated with 
perfect Andreev reflection, indicates the presence of a 
single localized Majorana zero-energy mode at the wire 
eilc J34jt38| p or a sufficiently high tunnel barrier, the con¬ 


ductance spectra will be peaked with this value at zero 
bias. While experimental results clearly show the devel¬ 
opment of such a peak upon tuning the system, at a finite 
magnetic field, into the predicted topological regime, the 
value of the zero-bias conductance peak is much less than 
the expected quantized value. The reasons for this dis¬ 
crepancy are addressed in RefsPMS anc [ alternative non¬ 
topological explanations have been advanced 41 ®! The 
lack of quant izati on of the experimental observations can 
be reconciled 5 with the Majorana theory by includ¬ 
ing the finite temperature and the finite length of the 
nanowire (thus allowing the Majorana modes from the 
two ends to overlap), but this physics is beyond the scope 
of our work where we restrict to zero temperature and a 
single normal metal-superconductor (NS) junction (as¬ 
suming the other Majorana mode to be far away from 
this junction). 

The difficulty in interpreting the tunneling experi¬ 
ments has prompted numerous theoretical studies on the 
co nduct ance of the nanowire device, using both numeri- 
ca j39iSH52 anc j analytical techniques®®. Although the 
latter works consider highly idealized models of the sys¬ 
tem, they are nevertheless valuable as they give clear 
insight into the parametric dependence of the trans¬ 
port physics as well as its dependence on various phys¬ 
ical properties of the experimental setup, which can 
then be applied to understand the more complicated nu¬ 
merical studies. An important question concerns the 
change in the conductance as the system is tuned from 
the topologically trivial to the nontrivial regimes (e.g., 
by tuning the applied magnetic field in semiconductor- 
superconductor hybrid structures). Remarkably, this as¬ 
pect of the physics has attracted relatively little atten¬ 
tion using these analytic method:!®. The purpose of this 
paper is to analytically address this aspect of Majorana 
physics in topological nanowire junctions. 

In this paper we examine the conductance spectra of 
one-dimensional NS junctions involving topological su¬ 
perconductors across their topological transition. We uti- 
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lize the Blonder-Tinkham-Klapwijk (BTK) formalism^, 
which is commonly employed to study junctions with un¬ 
conventional superconductonP^ES, to obtain analytic re¬ 
sults for the tunneling conductance of two models of a 
topological superconductor junction: a junction between 
a spinless normal metal and a p-wave superconductor, 
and a junction between a spinful normal metal and a 
spin-orbit-coupled s-wave superconductor in a magnetic 
field. The former is the simp lest m odel for tunneling into 
a topological superconductoJ^ 56 * 57 ( while the latter is a 
minimum modeP for the semiconductor nanowire device 
where experimental signatures for Majorana zero modes 
have been reported through the observation of zero-bias 
tunneling conductance peaks at the NS junction. We 
note that the spinless p-wave superconductor can be re¬ 
garded as an effective low-energy theory for the semicon¬ 
ductor nanowire, but this is inadequate for understanding 
the conductance spectroscopy of the device. Our analy¬ 
sis is analytical, and in particular we give explicit ex¬ 
pressions for the zero-bias tunneling conductance at zero 
temperature, which clearly shows an abrupt change at 
the topological transition. Specifically, we find that in 
the topological regime, the zero-temperature zero-bias 
conductance is quantized at a value of 2 e 2 /h indepen¬ 
dent of the barrier strength Z, but the detailed struc¬ 
ture (e.g., the width and the shape) of the quantized 
zero-bias conductance peak is controlled by the barrier 
transparency and the magnitude of spin-orbit coupling. 
Our BTK theory for the topological NS junction also 
shows that a finite barrier transparency could lead to the 
experimentally observed soft gap which is ubi quito us in 
semiconductor nanowire tunneling experiment^ 9 12 ~ 14 . 

The paper is organized as follows. In Sec. [TTJ we warm 
up by studying the conductance of a junction between a 
spinless normal metal and a spinless p -wave superconduc¬ 
tor across the topological transition. We then generalize 
the theory to consider the semiconductor nanowire device 
in Sec. Ill In particular, we obtain analytic results for 
the conductance spectra in the limits of a strong Zeeman 
field and strong spin-orbit coupling. Finally, the results 
are summarized in Sec. IIVI with a conclusion. 


II. JUNCTION WITH A SPINLESS p-WAVE 
SUPERCONDUCTOR 


the junction are generally unknown, and so Z is treated 
as a free fitting parameter. The Hamiltonian in each re¬ 
gion is written Hj(x) = \ f dx^(x)Hj(x)'i/j(x), where 
'kj(x) = {tlj^(x),il)j{x)) T are Nambu spinors and ipj(x) 
(ipj(x)) denotes the creation (annihilation) field oper¬ 
ator in region j = N (NM) and p (pSC). Assuming 
that the mass m is uniform throughout the system, the 
Bogoliubov-de Gennes (BdG) Hamiltonians are 

H n {x) = ( -h 2 d x /2m - p, N )r z , (la) 

K p {x) = {-h 2 dl/2m - h p )t z - iA p d x r x , (lb) 

where fix (p p ) is the chemical potential of the NM (pSC), 
A p > 0 is the p -wave pairing potential, and t p are the 
Pauli matrices acting on the particle-hole space. 

For notational simplicity, in the following we work 
with units such that H , /z^r, and 2m are all equal to 
unity. The energy spectra of the NM and pSG are 
then given by eN,±(k ) = ±(k 2 — 1) and ej L ±(k) = 
± \J ( k 2 — /ip) 2 + (A p k) 2 , respectively. In Fig. Ill we plot 
the spectrum of the pSC for different values of/Zp. Note 
that the spectrum becomes gapless at p p = 0 which 
marks the topological transition^ between the BEC-like 
strong pairing phase (p p < 0) and the BCS-like weak 
pairing phase (p, p >0). In the latter case, the positive en¬ 
ergy spectrum only develops the characteristic “double- 
well” BCS structure for p p > A 2 /2, with minimum value 

Ei = Apy//q, — A 2 /4 at k = ± Jp p — A 2 /2, and a local 
maximum value £2 = p p at k = 0. 



We start by considering a one-dimensional junction be¬ 
tween a spinless normal metal (NM) and a p-wave su¬ 
perconductor (pSC), which are located at x < 0 and 
x > 0, respectively. Their interface at x = 0 is mod¬ 
eled by a (5-function barrier of strength Z following the 
BTK prescription. The parameter Z controls the bar¬ 
rier transparency at the NS interface, and is the key 
parameter in the theory quantifying the tunneling con¬ 
ductance properties at the junction: a low (high) value 
of Z corresponds to a barrier with high (low) trans¬ 
parency at the NS interface. A microscopic evaluation 
of Z is typically difficult since the microscopic details of 


FIG. 1. (Color online) Typical energy spectra of the spinless 
pSC illustrating the nontopological (p, p = —0.01), transition 
(p, p = 0) and topological regimes (p, p = 0.01). In all curves 
we set A p = 0.05. 

We consider the scattering of an electron injected from 
the NM into the pSC with energy E. The incident 
electron can be normal reflected as an electron, An¬ 
dreev reflected as a hole, or transmitted into the pSC. 
The scattering wave function is $(a;) = ^n{x)Q(—x) + 
$p(x)0( x), where 0(x) is the Heaviside step function 
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and 


of the pSC wave function are solutions of the equation 


<M*) = (l) e “+(S) e -“» ( 2a ) 

$ p (x) = c ^ 7 -) elk ~ X + d ( 7 1 + ) eik+X ’ ( 2b ) 

where a and b are the Andreev and normal reflection 
amplitudes, respectively, c and d are the transmission 
coefficients into the pSC, and 


7± 


E T k± — p p 
A p k± 


( 3 ) 


Note that we approximate the wave vector of the elec¬ 
trons and holes in the NM by the Fermi momentum 
kp = y/2mprj/H, valid for E <C 1. The momenta k± 


E 2 = ( k 2 - ppf + (A p k) 2 . (4) 

Depending on the energy E of the incoming electron and 
the chemical potential p p , the pSC wave function can 
either be evanescent with complex solutions of Eq. 0 , or 
involve propagating states corresponding to real solutions 
of Eq. 0 with positive group velocity. We classify the 
different solutions in Table Q] 

The wave functions satisfy the continuity equation 
< J ) p( ;c )|x=o+ = < f > 7 v(a ; )U=o- and current conservation con¬ 
dition Jp<f> p (:r)| x=0+ - JAr4>Ar(a;)| x=0 - = - 2 iZr z $ A ,( 0 ) 
where the current operators are given by 

Jn — 2id x r z , Jp 2i3 x T z T A p t x . (5) 

Solving the boundary conditions, we derive the Andreev 
and normal reflection coefficients 


d{E) 

HE) 


A P (7+ - 7-) - 2(fc + - fc_) 

V E 

(2 - 2 iZ - k+ - k- - P)(7_ - 7+) + %(fc+ - fc-)(7+7- + 1) 

V E 


and 


( 6 a) 

( 6 b) 


respectively, where 


VL = 1 + (Z — ik-)(Z — ik+) - -E, and (7a) 

V E = D(y_ - 7 +) - ^(fc+ - fc_)(7+7_ + 1) - (k + - fc_)( 7 _ + 7 +). (7b) 


Pp 

E 

k -, fc+ 

Pp < A p /4 

0 < E < E/2 

E > E 2 

ki-, ki+ 
ki+, 

A p /4 < p p < A p /2 

0 E± 

Ei < E < E 2 

E > E 2 

kc~, kc+ 
ki-, ki + 
ki+, kp + 

Pp > A p /2 

0 < < Ei 

Ei < E < E 2 

E > E 2 

kc--, kc+ 
kR~, k R+ 
ki+, kp + 


TABLE I. Various solutions of Eq. Q for different val¬ 
ues of chemical potential p p and energy E, where E\ = 
A p \J p p — A 2 /4 and E 2 = \p P \- We denote propagating solu¬ 
tions by kn±, while evanescent solutions are given by ki± 
and kc±- These are given by kR± = ±[(p p — A p /2) ± 
VE 2 - El} 1 ' 2 , ki± = t[(AjS/2 - Up) ± VE 2 - E 2 ] 1 ' 2 and 
k c ± = ±[(p P - Ap/2) ± iy/E 2 - E 2 ] 1 ' 2 . 

Within the BTK formalisnP^the zero-temperature dif¬ 
ferential conductance is given by 

G p (E) = G 0 (l + \a(E)\ 2 -\b(E)\ 2 ), (8) 

where Go = e 2 /h is the normal state conductance for a 


quantum point contact. Although the general form of 
G p (E) is lengthy and unenlightening, relatively simple 
expressions can be found for the physically interesting 
case of zero bias, i.e., E = 0, which we provide in Table 
In] for the three different regimes of p p . In particular, we 
find that the zero-bias conductance abruptly jumps from 
G p (0) = 0 in the trivial regime (p p < 0) to G p (0) = 2 in 
the topological regime (p p > 0). The quantized conduc¬ 
tance is characteristic of the topological state, and can be 
interpreted as indicating perfect Andreev reflection [i.e., 
|a(0 )| 2 = 1, |6 (0 )| 2 = 0] at an interface supporting a Ma- 
jorana mod^SMEU it is therefore independent of the 
barrier strength Z and p-wave pairing potential A p . At 
the transition point (p p = 0) we find G p (0) < Go, with 
the exact value depending upon Z and A p . 

Characteristic plots of the conductance as a function 
of the energy are shown in Figs. [2] and [3] In general, 
the tunneling conductance G P (E ) decreases with the bar¬ 
rier strength Z, although in the topological regime the 
zero-bias conductance is unaffected by Z. Furthermore, 
it is interesting to note that in the topological regime, 
the width of the zero-bias peak decreases with Z but 
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f-Lp ■'C 0 

■ft 

■e 

II 

o 

/2p 0 

°(o) 

0 

2 Ap 

( Z + Ap/2) 2 + 1 + Ap 

—i 

6(0) 

—e ilp 

{Z + A p /2)“ + 1 itp 
( Z + Ap/2) 2 + 1 -f- Ap 

0 

G P (0) 

Go 

0 

[(Z + A p /2) 2 + l] 2 — A p 
[{Z + A p /2) 2 + 1 + A p ]2 

2 


TABLE II. Explicit expressions for the zero-bias Andreev 
reflection coefficient a(0), normal reflection coefficient 6(0), 
and differential conductance G p (0) for the spinless NM-pSC 
junction. The results are classified according to the three 
different regimes of fi p : the nontopological state (p p < 0), 
the topological phase transition point (/ i p = 0), and the 
topological state (p p > 0). The quantity <p is defined by 
sin ip = 2(Z + yj A p /4 — [j, p )/\(Z + \J A p /4 — fi p ) 2 + 1]. 

fi p = — 0.01 fj.p = 0 lip = 0.01 



E E E 

A p = 0.07 . A p = 0.7 - A p = 2 - A p = 20 


FIG. 2. (Color online) Variation of the tunneling conductance 
G P (E) with pairing potential A p and chemical potential p, p 
for the spinless NM-pSC junction. The values of the pairing 
potential A p are given in units of /xjv/fcp, while the chemi¬ 
cal potential p, p and energy E are expressed in units of /in- 
We show typical results for the nontopological (/r p < 0, left 
column), transition (p p = 0, middle column), and topologi¬ 
cal (p p > 0, right column) regimes, and for barrier strength 
Z = 0 (top row) and Z = 1 (bottom row). 


shows a nonmonotonic dependence with A p ; the width 
first increases as A p increases, however, beyond a cer- 



E E E 

Z = 0 - Z = 1 . Z = 2 


FIG. 3. (Color online) Variation of the tunneling conductance 
G P (E) with the barrier strength Z and the pairing poten¬ 
tial A p for the spinless NM-pSC junction in the topological 
regime. The values of the pairing potential A p are given in 
units oi (in/ kF while the chemical potential /.t p and energy E 
are expressed in units of /in- Note that the zero-bias conduc¬ 
tance is constrained to be 2Go by the topological condition. 


tain value of A p , the width decreases with A p . For 
Up < Ap/2, a singularity appears in the G P (E) curve 
at the gap edge E 2 = H P \ ■ On the other hand, as shown 
in Fig. |3])a), two singularities are visible in the conduc¬ 
tance for fi p > Ap/2, corresponding to the edge of the 

gap at Ei = A p yjp p — A p /4 and the local maximum in 
the spectrum at E 2 = /i p . 


III. JUNCTION WITH A 
SPIN-ORBIT-COUPLED NANOWIRE 

In its topological phase, the low-energy sector of the 
spin-orbit-coupled nanowire proposal is formally equiv- 
alent to the spinless p- wave superconductor studied 
abovd^ 7 ^. In order to obtain the conductance spectrum 
and its variation across the topological transition, how¬ 
ever, we must examine the full model including spin-orbit 
coupling and Zeeman splitting. In this section we there¬ 
fore consider a one-dimensional junction between a spin- 
split spin-orbit-coupled superconducting wire (SOCSW) 
and a spinful normal metal (NM), which occupy the re¬ 
gions x > 0 and x < 0, respectively. Similar to Sec. [TTJ 
we model their interface at x = 0 by a 5-potential bar¬ 
rier of strength Z. The Hamiltonian in each region is 

written Hj(x) = \ f dx^j(x)'Hj^j(x), where \&j(ic) = 

and -ipj a (x) [ipja(x)] is 
the creation (annihilation) field operator of an electron 
with spin er in region j = N (NM) or S (SOCSW). Using 
the same unit convention as in the previous section, we 
write the BdG Hamiltonians of the NM and SOCSW as 

Un = (-8% - 1 )r z , (9a) 

Us = -d 2 x T z - iad x r z a z + V z cr x + A 0 t x , (9b) 
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momentum, k 




momentum, k 


FIG. 4. Energy spectrum of the SOCSW for different values 
of Zeeman potentials: (a) Vz = 0 (nontopological), (b) Vz = 
0.01 (nontopological), (c) Vz = 0.02 (transition), and (d) 
Vz = 0.03 (topological). In all plots, we set a = 0.3 and 
Ao = 0.02. 


where cr M (r p ) are the Pauli matrices in spin (particle- 
hole) space, a is the strength of spin-orbit coupling 
(SOC), Vz is the Zeeman field, and A 0 > 0 is the 
proximity-induced s-wave pairing potential which is as¬ 
sumed to be real. We set the chemical potential of the 
SOCSW and Zeeman coupling in the lead to be zero, and 
take uniform electron masses throughout the system. 

The positive branches of the BdG spectrum of the 
SOCSW are given by 


E ± =( k 4 +a 2 k 2 + A 2 0 + V z 


± 2\Jk 4 (a 2 k 2 + Hj) + AgPj j ' . (10) 


As shown in Fig. [4j the energy spectrum is gapped ex¬ 
cept for V z = A 0 . This value of V z marks the topolog¬ 
ical quantum phase transition between the topologically 
trivial (V z < A 0 ) and nontrivial phases ( V z > Ao(F" 7 '. 
Although Eq. (10) can be analytically solved for the mo¬ 


menta corresponding to a given energy E, the general 
expression is unwieldy. In what follows, therefore, we 
will instead work in the limits of a strong Zeeman field 
and strong SOC, where more compact results can be ob¬ 
tained. 


A. Strong Zeeman splitting 

In the limit of strong Zeeman splitting (V z a, Ao), 
the quasiparticle excitation spectrum of the SOCSW is 
split into two spin bands as shown in Fig. [5ja) . In the 
normal state (Ao = 0) the spectrum is approximately 
given by e±(k) ~ k 2 ±V z . The system is then a half¬ 
metal, with only one spin-polarized band [e_ (k)] occu¬ 
pied. Projecting the full Hamiltonian into this band gives 



FIG. 5. Energy spectrum of the SOCSW in the limits of (a) 
strong Zeeman field and (b) strong SOC. For clarity, only the 
positive energy branches of the spectrum are shown. In panel 
(b), the spectrum about the minima at k = 0 constitute the 
“interior” branches, while the spectrum about the minima 
at k = ±a are the “exterior” branches. Note the different 
effective gaps for these branches, and the states contributing 
to the slowly varying left- and right-moving fields, L a (x) and 
R a {x), respectively. 


the effective HamiltoniarP ®III 

H s ( k ) = {e-(fc)V’^_(fc)^s-(fc) 

}■ 

( 11 ) 

where A_(fc) ss akA$/V z is a p-wave pairing potential 
and annihilation (creation) field opera¬ 

tor for e_(fc) band. The projected Hamiltonian is equiv¬ 
alent to the spinless pSC Hamiltonian H p {k ) [Eq. 0|. 
with the identifications p p = V z and A p = aAo /V Z - If 
the Zeeman field is applied on both sides of the junc¬ 
tion such that the NM is also fully spin polarized, then 
the low-energy sector is identical to the spinless NM-pSC 
junction studied in Sec. [TTJ and the results obtained above 
for the differential conductance directly apply. 


A_(/c)-i/ig_(fc)^g_(— k) + h.c. 


B. Strong spin-orbit coupling 


In the case of strong SOC (a V z , A 0 ), the BdG spec¬ 
trum of the SOCSW has the characteristic form shown 
in Fig. [5jb). In particular, we note that both the + 
and — spectra [Eq. (10)] have minima at k = 0 (the so- 
called interior branches), while the — spectrum also has 
minima at k = ±a (the exterior branches). For small 
energies E < Aq, V z . we can linearize the Hamiltonian 
about these minima by introducing the ansatz for the 
field operatord^Sfil 


V’stM ~ + h( x ) e iax , (12a) 

ipsi(x) ~ L±(x) + Ri(x)e iax , (12b) 

where R a (x) and L a (x) represent slowly-varying right- 
and left-moving fields, respectively; see Fig. [5Kb) . In¬ 
serting this ansatz into the Hamiltonian [Eq. (|9b|)] and 
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neglecting all “fast oscillating” terms (involving terms 
with phase factors e ±lax ), we obtain effective Hamilto¬ 
nians valid for the states near the interior and exterior 
branches. Specifically, we write 


and 


vl + u 


2 

v 


1, for E > |A„|, 

E/ |A„|, for 0<E< |A„|, 


(18) 


H s ] = 1 
s 2 


I (13) 


where l = e, i denotes the exterior and interior branches, 
respectively, and the BclG Hamiltonians are written 


'H ( g ) = - iar z a z d x + A 0 t x , 

= -ioLT z o z d x + V z a x + A 0 t x . 


(14a) 

(14b) 


The spinors for the interior and exterior branches 
are defined in terms of the slowly-varying field as 

$5 e) ( x ) = R l( x )’ R \( x )i ~ L \( X )) T and ^(z) = 

(i? T (x), L±(x), L[(x), -i?|(x)) T . 

We consider an electron with energy E and spin a in¬ 
jected into the SOCSW from the NM. The wave function 
in the NM is given by 





dWaOz;) = 


( 15 ) 

where 5 a <j> is the Kronecker symbol. The coefficients 
a aa i and b a<J i are the amplitudes for Andreev and nor¬ 
mal reflection, respectively. Note that due to the SOC in 
the SOCSW, both spin-flip and spin-preserving reflection 
processes are allowed. The wave function in the SOCSW 
is a superposition of solutions on the exterior and interior 
branches 


®s<r(x) = 


/ -u_ \ 


( u +\ 

I sgn(A_)v_ ' 

e ik - x +c^ 2 | 

v+ | 

— sgn(A_)r>_ i 

r + J 

V u- ) 


\U+J 


+ c (e) 

' C <71 


f v A 

I 

0 

e i{k^- a)x + c (e) 

u 0 

\o) 

\ 


e ik + 3 




(16) 


where the first line on the right-hand side gives contri¬ 
butions from the interior branches, while the second line 
originates from the exterior branches. Note that the coef¬ 
ficients 2 ) and 2 , are the transmission coefficients 
into the SOCSW. The elements of the wave function are 
given by 


u 


2 

V 


(E+^EP-Al) /2 E, 

[E + i^Al-E*) /2|A„|, 


for E > |A„|, 
for 0 < E < |A„|, 
(17) 


where v = ±, 0, and A± = A 0 ± Vz- The wave vectors 
appearing in Eq. (16) are E 2 — Aj_/a for the in¬ 

terior branches, and = -JE 2 — Aq/q for the exterior 
branches. 

The wave functions satisfy the continuity and current 
conservation boundary conditions 


$s<7(z)U=o+ = $JVa(a;)U=o- , (19a) 

>/s < I ) Sct(^)U=o+ — Jn®n<t(x)\ x =o- = — 2iZT z $ Na (0), 

(19b) 


where the current operators are given by 

Jn = 2 i() :!: t z , Js = 2 id x T z -f- cxt z ct z . (20) 

In the limit of strong SOC (a>Vz, A 0 ), we ignore terms 

(i) (e) 

proportional to k_ , k\ , A?q -C 1 in the current conser¬ 
vation equation. Expressions for the Andreev and normal 
reflection coefficients found from solving these equations 
are given in Appendix [A] 

The zero-temperature differential tunneling conduc¬ 
tance Gs{E) is obtained from the BTK formula 


Gs{E) = Go ( 2 + Y, {ME)\ 2 - KdE)\ 2 } 
<7,€=t,4 

( 21 ) 

Although the general expression is complicated, compact 
forms for the reflection coefficients and the conductance 


at zero bias are presented in Table III As in the spinless 


NM-pSC junction studied above, the zero-bias conduc¬ 
tance Gs(0) is discontinuous across the topological phase 
transition. In the topological regime (Vz > A 0 ) the 
zero-bias conductance takes the quantized value Gs(0) = 
2. This implies that the Andreev reflection coefficients 
in Eq. © exactly cancel the normal reflection coeffi¬ 
cients; moreover, from Table |III| it can be verified that 

E<t,€ KfK°)l 2 = E<t,« l^(°)r == L This can be un¬ 
derstood in terms of the existence of a single Majorana 
mode at the interface which couples to one of the two 
channels in the normal regioipSHSI. While there is perfect 
Andreev reflection in this channel, in the other channel 
we have perfect normal reflection. In the nontopological 
regime, on the other hand, Gg(0) takes on nonuniversal 
values and is dependent upon Z and a. In particular, 
the zero-bias conductance in the nontopological phase 
can strongly exceed the quantized value in the topologi¬ 
cal state: for the gapped nontopological state ( V z < A 0 ) 
and at the topological transition point (Vz = Aq), we 
find the maximum values Gs(0) = 4 and Gg(0) = 3, 
respectively, which are realized for Z = 0 and a = 2. 

We plot the calculated conductance as a function of 
energy in Figs. [6] and [7J In the former we show examples 
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2 


TABLE III. Zero-bias values of the Andreev reflection coefficients a atT /{ 0), normal reflection coefficients b aa /(0), and differential 
conductance Gs (0) in the strong SOC limit of the NM-SOCSW junction. The three columns give the values in the nontopological 
(Vz < Ao), transition {Vz = Ao) and topological (Vz > Ao) regimes. The terms D 1,2 are given by D\ = 2[1 + Z 2 + {a/2) 2 ] 
and D 2 = Z 2 + (1 + a/2) 2 . 


of the conductance spectra in the nontopological, tran¬ 
sition, and topological regimes, while the latter explores 
more fully the variation of the conductance spectra in 
the topological state away from zero bias. The conduc¬ 
tance spectra show a much more complicated structure 
than those in the spinless NM-pSC junction, reflecting 
the presence of three distinct gaps (A + , |A_|, A 0 ) in the 
strong SOC limit of the SOCSW. Indeed, at the energy 
corresponding to each gap we observe a nonanalyticity 
in the conductance spectrum. Although there is consid¬ 
erable variation in the conductance spectrum as a func¬ 
tion of energy, a number of trends can be discerned: in¬ 
creasing Z tends to suppress the conductance, the energy 
variation of the conductance is nonmonotonic in general 
with cusplike structures at specific energies, and the en¬ 
ergy variation of the conductance is stronger near zero 
energy for larger values of Z. While the conductance at 
first tends to be enhanced by increasing the SOC, the 
conductance eventually goes through a maximum before 
monotonically decreasing. Similarly, the SOC increases 
the width of the zero-bias peak in the topological regime, 
but beyond a certain SOC strength it decreases again. 
The basic finding is that, other than the universal quan¬ 
tized Majorana peak at zero energy, the tunneling con¬ 
ductance shows interesting and nontrivial dependence on 
Z and E in the topological phase. In particular, au in¬ 
teresting conclusion of our theory is that the zero-bias 


conductance could be quantized in the topological phase 
for small values of Z without developing a peak in the 
tunneling conductance at all. 

Note that the above discussion holds true also for the 
case where the Zeeman coupling in the lead or the chemi¬ 
cal potential /rg of the SOCSW are nonzero. For the case 
where |A_| < Aq, the zero-bias peak formed in the topo¬ 
logical regime is within an energy range of | A_|. Since 

the topological gap |A_| = |^/Aq + — Vz\ decreases 

with the absolute value of the chemical potential |//g|, 
the width of the zero-bias peak decreases with |/ig|. 


IV. SUMMARY 

Using the BTK formalism we have analytically stud¬ 
ied the zero-temperature tunneling conductance spectra 
of NS junctions involving topological superconductors. 
Finite temperature effects within this formalism simply 
lead to thermal broadening of the zero-temperature con¬ 
ductance and can be included in the theory numerically 
by introducing an integration over the Fermi function®'. 
As in the BTK papef®^ the finite-voltage conductances 
are found to depend on the strength of the barrier at the 
interface, which is parameterized by the dimensionless 
parameter Z. Specifically, we have examined a spinless 
NM-pSC junction and a spinful NM-SOCSW junction, 
































Vz = 0.5A 0 


Vz = A 0 Vz = 1-5A 0 



FIG. 6. (Color online) Variation of the tunneling conduc¬ 
tance Gs(E) with SOC strength a and Zeeman field Vz in the 
strong SOC limit of the NM-SOCSW junction. We present 
typical results for the nontopologic.al (Vz < Ag, left col¬ 
umn), transition (Vz = 0, middle column), and topological 
(Vz > Ao, right column) regimes, and for barrier strength 
Z = 0 (top row) and Z = 2 (bottom row). In all plots we set 
Ao = 0.001. The values of Ao and Vz are given in units of 
Hi v, while the values of a are expressed in units of Hi v/fep- 


paying particular attention to the change in the zero-bias 
conductance across the topological phase transition. We 
explicitly demonstrate that the zero-temperature zero- 
bias conductance is quantized at a value of 2 e 2 /h in the 
topological regime, in agreement with effective models of 
these systems based on a single Majorana mode coupled 
to a normal channel. Despite this quantization at zero 
voltage, the zero-bias conductance only develops a peak 
(or a local maximum) as a function of voltage for barriers 
with sufficiently large Z parameter, or for small and large 
SOC strength. These parameters also control the width 
of this peak. In the nontopological regime, on the other 
hand, the conductance takes nonuniversal values depend¬ 
ing upon the details of the system. In both cases the 
conductance spectrum away from zero bias shows con¬ 
siderable variation with the details of the junction. Our 
calculated BTK conductance also shows that the con¬ 
ductance is finite inside the superconducting gap region 
because of the finite barrier transparency, providing a 
possible mechanism for the observed “soft gap” feature 
in the experimental studiepb-dfH] This effect is qualita¬ 
tively similar to the “inverse proximity effect” at the NS 



E/Ag E/ Ao E/ A 0 

-Z = 0 - Z = 1 . Z = 2 


E/Ag E/Ag E/Ag 

-Z = 0 - Z = 1 . Z = 2 


FIG. 7. (Color online) Variation of the tunneling conductance 
Gs(E) with interface barrier strength Z and spin-orbit cou¬ 
pling strength a for the strong SOC limit of the NM-SOCSW 
junction in the topological regime. We set Aq = 0.001 and 
Vz = 1.5Ao. The values of Ao and Vz are given in units of 
Hn, while the values of a are expressed in units of hn/^f- 
Note that in all cases the zero-bias conductance is equal to 
2Go, consistent with the topological state. 


interface arising from the finite barrier at the interface as 
discussed in the recent literature 50 !, although other pos¬ 
sible physical mechanisms for the soft gap behavior have 
also been proposed^. We mention finally that our theory 
is for a single NS junction which effectively assumes the 
existence of only a single Majorana mode at the NS in¬ 
terface (with the other Majorana bein g loca ted infinitely 
far away) and thus Majorana splittin^ 45 * 69 ^ 1 due to the 
wave function overlap between two Majorana modes is 
not germane to our theory (but can be included if neces¬ 
sary in a future generalization). 

Finally, we emphasize that one of the most salient fea¬ 
tures of our theoretical work is that it is completely an¬ 
alytical within a continuum model in contrast to most 
theoretical works on Majorana nanowires which focus 
on numerical simulations within a tight-binding lat¬ 
tice model. All the microscopic details of the complex 
normal-superconductor tunneling process are simply sub¬ 
sumed in a single phenomenological parameter Z (“the 
interface barrier strength” of the BTK formalism) allow¬ 
ing our theory a great deal of flexibility for actual model¬ 
ing of the Majorana nanowire experimental results, since 
the realistic microscopic details of the NS interface are 
rarely known in the actual experimental nanowire setups. 
It is gratifying that our analytical model captures the es¬ 
sential features of the Majorana nanowire experiments 
through our finding of the Majorana zero-bias conduc¬ 
tance quantization and soft gap feature, with the inter¬ 
esting prediction that for strong metallic junction (i.e., 
for very low interface barrier or a small value of Z) the 
Majorana zero-bias conductance may not necessarily be a 
peak in the tunneling conductance although it would still 
be quantized since the Majorana zero mode necessarily 
implies perfect Andreev reflection. 
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Appendix A: Andreev and normal reflection coefficients for the NM-SOCSW junction 


Solving the boundary equations, we obtain the Andreev (normal) reflection coefficients a CTCT ' {b a , T >) as 

au 0 v 0 [u-v+ — sgn(A_)u + u_][l + (Z + if) 1 2 3 4 5 6 * 8 9 ] 

=-WiM2j-’ 

U E U E 

V+ / \ \ v — 

— TT +sgn(A_)- 7 . 


a^i(E) = au 0 

+ j). 

au 0 v o [u-V + - sgn(A_)u + u_] [l + (Z - if) 2 ] 


au( E ) = -- 


x4 1} x4 2) 


b n (E) = b u {E) 

= ~ _ ( * + Z)2+ (f) 

| sgn(A_)ugU_u + [Z 2 + (f - l) 2 ] - u 0 v 0 [u-v + + sgn(A_)u + u_][l + Z 2 + (f ) 2 ] + ulu-U + [Z 2 + (1 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 


f) 2 ] 


V E )V{ E 


, au 2 0 {u-v+ - sgn(A_)u+u_](l - iZ + f ) 2 

bn( E ) = - 


v( e v( e 


bit(E) = 


][m_u + - sgn(A_)u + u_](—1 + iZ + f ) 2 


v{ e v( e 


(A6) 

(A7) 


where = uqu + [Z 2 + (a/2 + l) 2 ] — VqV + [Z 2 + {a/2 — l) 2 ] and = uqU- [Z 2 + {a/2 + l) 2 ] — sgn(A_)t>of 
(a/2-1) 2 ]. 
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